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1 Introduction 

The annihilation cross section is one of the key ingredients of computing the thermal relic 
density of the dark matter, as well as the indirect detection signals. Due to the fact, that 
until recently the observational data for the relic density and cosmic rays spectra were 
far from being accurate, most of the work in the literature uses only its tree level value. 
In recent years however, some work has been done in order to go beyond and include 
one-loop corrections [1-3], as well as Bremsstrahlung processes [4-15]. It has been found 
that including higher order processes can lead to rather large corrections and also may 
significantly alter the cosmic rays spectra. 

Another effect studied recently in this context is the so-called Sommerfeld enhancement 
[16], that is the modification of the wave function of the incoming non-relativistic particles 
due to their mutual interaction (which is non-perturbatively treated). It has been shown to 
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be a very important modification in a number of models with new dark forces (see e.g. [17] 
for a recent review), and even in more standard cases like in the Minimal Supersymmetric 
Standard Model (MSSM) [18-23]. 

In this work we describe how to incorporate these corrections simultaneously and give 
the full one-loop and Sommerfeld corrected results for the case of the dark matter being a 
fermion living in the adjoint representation of SU{2)wr. This is the case of the pure Wino 
neutralino, but it is also interesting case per se (see e.g. the Minimal Dark Matter model 
[24]). It is also a good starting point for possible extensions because its relative simplicity 
makes more clear the description of the various effects. Another, more phenomenological 
reason to study this case is that it is precisely the one in which the Sommerfeld effect in 
the MSSM is the most important. 

We start the discussion from presenting the model and reviewing the Sommerfeld 
effect. In section 2 we compute the radiative corrections and present the results. In section 
3 we compute the full cross-section (including the Sommerfeld enhancement). Finally we 
very briefly mention the possible phenomenological applications and give our conclusions 
in section 4. 

1.1 The model 

We consider a model in which a Majorana fermion plays the role of the dark matter. We 
call it and assume that it belongs to the adjoint representation of the SU{2) subgroup 
of the electroweak SU{2) x U{1). The other two members of the triplet can be combined 
to be described as a charged Dirac fermion and its anti- fermion which we call x^- We 
assume that this triplet of fermions is massive due to an explicit mass term, that is present 
independently of the Higgs mechanism that might give mass to the weak vector bosons, 
and in fact we assume that these fermions do not interact with the Higgs field (if any). 
Within this setup, interacts only with the charged weak-interaction vector bosons 
and its charged partners only interact with the W^, and with the Z and 7. 

We will be interested in the the mass m of the Wino up to a few TeV. For the computa- 
tion of the radiative corrections we also assume that the charged fermions of the multiplet 
have a mass which is higher of a negligible amount with respect to the TeV scale, however 
we will include this difference in the Sommerfeld effect.^ 

In our model, at the tree level there is only one possible annihilation channel:^ 

^0^0 ^W+W-. (1.1) 

However, at a higher order it is possible that the x^X^ paii' becomes a (real or virtual) 
X~^X~ pair which subsequently annihilates: 

X^X'^W+W- or x^X^ ^ ZZ,Z^,jj. (1.2) 

^The mass difference comes from radiative corrections and is of the order of 5m = 0.2 GeV [25]. 

^In the pure Wino scenario in the MSSM there are additional annihilation channels. However, in the 
case in which we are most interested in, i.e. x° having a mass in the TeV range, this channel is by the far 
dominant one. 
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This sequential process is formally of higher order, but it can be enhanced by the two- 
channel version of the Sommerfeld effect (see section 1.2). In this case it can be effectively 
of the same as the tree level process and has to be included. 

Moreover, since we consider the radiative corrections which provide an order 0{g'^) 
correction to the amplitude and corresponding to an order 0{g^) term to the cross-section, 
we have also to include the annihilation in three final particles, i.e. 

// ^ W^W'Z, W^W~-i (1.3) 

and, by Sommerfeld effect, also 

X°X° ^ X+X" ^ W^W-Z, W+W-^ . (1.4) 

We are interested in the case when annihilating particles are non-relativistic, therefore, 
we can take ctq to be the nominal cross-section for the annihilation at rest within a negligible 
relative error 0{v'^), that is the relative variation of the Mandelstam variables s and t 
averaged over the angles. In this case the two incoming neutralinos, being Majorana 
fermions, form an s-wave spin-singlet. This is a very good approximation for dark matter 
particles in the halo today, and it allows for a great simplifications of the computations. 
Firstly, because in this case the initial pair (being Majorana fermions) have to be in a 
s-wave spin singlet. Secondly, the kinematics simplifies, since the annihilation becomes like 
a decay of a particle with the mass 2m. 

This approximation however sets limits on the usage of the results for the relic density 
calculations. Although at freeze-out the dark matter particles are still non-relativistic, 
their velocity is about v ~ 0.3. On the other hand, in order to get accurate results in this 
case one needs not only to generalize this computations, but also include the p-wave, which 
is beyond the scope of this work.'^ 

1.2 The Sommerfeld enhancement 

The Sommerfeld enhancement (for a recent review see [22] and references therein) of a 
process is usually stated to give the full cross-section of the process in terms of the multi- 
plication of two factors: 

a = S{v)aQ (1.5) 

where S{v) is the velocity-dependent Sommerfeld factor, which represents the possible 
increase (decrease) of the fiux of the incoming particles, due to their mutual attraction 
(repulsion) and ctq is the nominal cross-section. Typically as uo one uses the tree level 
value, but if one wants to incorporate loop corrections, it can be also computed at higher 
order in perturbation theory.^ In our work we are going to compute ctq up the order 
0{g^). Whereas we compute the nominal annihilation at rest, we still keep into account 

^For some results including one-loop corrections to relic density computations, however without the fully 
treated non-perturbative Sommerfeld effect see refs. [2, 3]. 

''This is true as long as the annihilation process is short distance one, so that the long distance Sommerfeld 
effect is decoupled and can be treated separately. 
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the velocity in the Sommerfeld enhancement, since in this process the momentum transfer 
can depend substantially on it. 

The Sommerfeld enhancement factors are computed by a formalism using a system of 
non-relativistic Schrodinger equations. In the case at hand, there are two possible channels 
through which the annihilation can take place: (x^^X^) ^iid The x^ P^-ii" can 

only become a (real or virtual) x^ P^^^ by exchange of W^, whereas the x^ P&i^ can 
self-interact by exchanging Z or photon and also come back to the x^ P^ir by exchange of 

This dynamics can be represented in terms of Feynman diagrams (see ref. [26]). For 
small velocities, it is well described by the non-relativistic evaluation of ladder diagrams, 
where the steps of the ladder are either W^, Z or 7 and the lateral bars are x^ or x^- 
In the non-relativistic approximation (neglecting spin-orbit effects) the total spin of the 
two-body system is conserved, therefore both the pairs x^X^ X~^X~ in ^ s-wave 
spin-singlet. Summing the ladder diagrams is equivalent to solving the coupled Schrodinger 
equations for the two body (reduced) wave functions 

Since in computing one-loop radiative corrections there are also diagrams containing 
the exchange of W in the form of a single ladder, one has in this case to subtract the 
non-relativistic part which is already included in the solution of the Schrodinger equations. 

2 

For large m, this non-relativistic part can be quite large C>( j:^-^^). The solution of the 
Schrodinger equations represents a non-perturbative re-summation of those large terms 
and this is one reason for including the Sommerfeld effect. 

We call ip^{x) and f^{x) the s-wave reduced wave functions for the xPx^ and x'^X~ P^-ii' 
respectively {x = rp, p = mv, v being the x^ velocity). In the low velocity approximation, 
the spin-singlet x^X^ state is described by 



a| |(6| |) being particle (anti-particle) creation operators at rest, for a given spin-projection. 
Note the identities, for any Dirac matrix M: 



/(x)ala!|0) 




and the spin-singlet x^X state by 
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The coupled equations are {£ = p^/m is the energy of the incoming pair and 5m = m±—m): 



a2(^0(x)V(x) + ^^^v^{x) = 0, (1.10) 



25m\ ^ 
hp (x) 




W X 
m V 



+ 4 ^^(x) + V2 if'^ix) 



V X 



(1.11) 



Here q = "^j^'" , = cos0^ and = sinOyj, with 9y^ being the Weinberg angle and 
a = 5'^/(47r) where g is the SU{2)l coupling. We take the values at the electroweak scale 
computed in the MS scheme [27]: = 0.23116 and g = 0.65169. One has to require 
the boundary conditions such that the (/J^ is incoming and outgoing, whereas tp^ is only 
outgoing (if 26m < £) or exponentially decaying (if 25m > £). 

As we explain below, in our case velocity is too small to allow for on-shell ip^, hence 
the relevant case is only the latter one and the boundary conditions are 



lim (i/(a;) - d^^°{x)) = -e"*^, lim ip^{x) = 0. (1.12) 
With this normalization, we call the amplitude Sommerfeld factors: 

So = dx^p^{x)\x=o, s± = dx^p^ix)\x=o, (1.13) 

then the (Sommerfeld enhanced) amplitudes of the annihilation processes for any Standard 
Model (SM) final state are: 

In the literature one often finds the one channel version of the Sommerfeld. In this case 
one would have a single wave function ip{x) and the Sommerfeld enhanced the cross-section 
would be given by eq. (1.5) with 

S'(-y) = |sp, s = d:cV{x)\;,=o. (1.15) 

However, if there are more channels, one has to use the formula for the enhanced amplitude 
eq. (1.14). Note that by taking the modulus square of eq. (1-14) there is also a cross term, 
which was neglected in the previous works on the Sommerfeld effect. 

The result for the Sommerfeld factors depends strongly on the mass-splitting between 
X° and x^'i in fact, also the computation is somewhat different if the total energy of the 
pair is greater or smaller than twice the mass. However, taking the mass splitting to be 
of the order of 0.2 GeV, and the velocity of the order of present day dark matter velocity 
V ~ 10~^, even for m being a few TeV, the production of real x^ from the x^ pair is not 
allowed. In this case (i.e. far below the x~^X~ threshold) so,± nearly does not depend on 
the velocity [22, 28]. 
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It is worth noting, that for small velocities the Sommerfeld effect at one loop level (i.e. 
not summed over all orders), gives: 

2 

so = l, s± = V2^^. (1.16) 
47r mw 

This result is recovered from the full solution of the Schrodinger equations when one takes 

2 

the limit of small v-^ and small -f — —. However, since for large enough m this does not 

m-w 47r mw ^ 

hold, the full numerical computation is needed. 
2 Electroweak corrections 

In this section we will present the computation and the results for the radiative corrections 
to the annihilation amplitudes. In order to compute the cross section up to the order 0{g^) 
we need to include both the one-loop corrections and the real production. In particular, 
the Bremsstrahlung of a soft /collinear photons is crucial for the cancellation of the IR 
divergences, as we will discuss in detail in section 2.2.4. 

We will organize the results in the following way. Firstly we discuss the issue of at 
what scale should we take the couplings. Then we consider the radiative corrections to the 
amplitudes A^o^o^\y+]^- , A^+y^-^iy+^y- , A^+^-^2Z,z^,^"f 

2.1 The scale of the coupling 

It is very important to understand at what energy scale the SU{2) coupling g and of the 
Weinberg angle 9^ should be taken. We argue, that both in the radiative corrections and 
the Sommerfeld effect computations we should take them at the electroweak (EW) scale. 

In the computation of the Sommerfeld factors so,± one could wonder whether one 
should take g at the scale of the neutralino mass m, since this sets the energy scale of the 
process. However, what matters in the Sommerfeld enhancement computation is in fact 
the scale of the momentum transfer between the incoming particles and this can be at most 
of the order of the vector boson mass. 

The detailed computation, both analytical and numerical, shows [29] that the radiative 
corrections to the vertices XX^^ with x on-shell as appropriate in the non-relativistic case, 
at zero momentum transfer exactly compensate the effect of the x wave-function renormal- 
ization. Therefore, there is no dependence on m once taking the renormalization scale to 
be mw, and this compensation persists quite effectively up to momentum transfer 0{mw)- 
Therefore, since what remains to be considered is the W wave-function renormalization, 
which is already included in the definition of the coupling at the scale mw, there are no 
appreciable corrections at all, if we take g at the scale mw- 

In fact, note that the use of the running coupling constant is appropriate for the 
processes which depend significantly on a single large scale. For instance it would be 
appropriate to take the coupling at the scale m for processes in which the momentum 
transfer to x is also of the order of m. 

In the case of the radiative corrections to the annihilation amplitude, there is not a 
precise compensation of the radiative correction of the vertices XX^ with the x wave- 
function renormalization because the internal x lines are off-shell, and therefore we take 
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into account these loops that give a further radiative correction not included into taking g 
at the scale mw- 

Note that, by computing the Feynman diagrams giving the vertex corrections and 
the wave function renormalization and fixing the renormalization at the EW scale, we are 
evaluating perturbatively how the coupling "runs" from its EW value. 

This is like expanding at the one-loop order the formula for the running coupling 
constant, except that we do not have to include the W wave-function renormalization, 
because it only depends on the square Vl^-four-momentum which is equal to and 
therefore it is already inside the definition of the coupling at the scale mw- 

Let us also recall that the standard use of the renormalization group techniques holds 
in the "deep euclidean region" in which the external lines are quite off-shell. In our case 
instead, the external particles are on-shell and therefore there occur not only the large log's 
related to the UV divergences but also large log's due to IR effects. As we will discuss in 
the following, we do not attempt a re-summation of the large log's of various origin, and 
this is another reason why we do not attempt to use a kind of non-perturbative formula for 
the running coupling, suitably modified to take off the W-wave function renormalization, 
which would correspond to some partial re-summation of one subset only. 

2.2 The radiative corrections to x^X^ ~^ W~^W~ 

We start the discussion from the one-loop corrections to x^X^ annihilation. Firstly we will 
discuss the method of doing the computations and the in section 2.2.3 we will give the 
results. The way we present them is in terms of the correction to the tree level amplitude: 



A = A 



tree 



1 + 



9' 



(47r): 



:Ci{ni] 



(2.1) 



where Ci{m) are the coefficients corresponding to the diagram i. 
2.2.1 The UV divergent diagrams 

The UV divergent one-loop diagrams come from the vertex corrections and the fermion 
wave- function renormalization, as presented on figure 1. 



w+ 




Figure 1: The UV divergent diagrams for x^x" ~^ W^W process. The vertex corrections 
(diagrams 1 and 2) and the fermion wave-function renormalization (both diagrams are included in 
3). 

For all these diagrams, we have done the computations using full analytical expressions 
with Feynman parameters and integrated analytically (using Mathematica) over the first 
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one and numerically over the second one.^ We took the Feynman-'t Hooft gauge for the 
W propagator, which simplifies the computations, noting that x is not coupled to the 
Higgs bosons and that there is no vertex with two Ws and one neutral unphysical Higgs. 
We used dimensional regularization and dropped the terms 0{l/e) because they are taken 
into account in the renormalization at the scale miy. In fact, the loop corrections to the 
coupling g evaluated at the mw scale do contain the same 0(l/e) terms, which therefore 
are part of the definition of the coupling at that scale. We also didn't include the W 
wave-function renormalization of the final VF's for the same reason.® 



2.2.2 The UV finite diagrams and the IR divergence 

Besides the loops giving the radiative correction of the vertices and the x wave-function 
renormalization, there are two other loops, which are not UV divergent.^ 






X 




Vvw 








x° x" 





4) 5) 
Figure 2: The UV finite diagrams for x°X° ~^ W~^W~ process. 



Diagram 4 represents a process in which the incoming x P^^^ goes to a virtual x 
pair (which then annihilates in W^) by exchange (see figure 2, diagram 4). The 
contribution of this loop is very large when m/mw is large. In fact, it is recognized that 
it contains the first order contribution to s±.^ This is seen because, as shown in ref. [26], 
the Sommerfeld effect comes by summing the non-relativistic part of the ladder diagrams, 
and this diagram is precisely the first of the series. 

More in detail, the statement for the Sommerfeld enhanced amplitude A = sAq comes 
from taking the amplitude as the non-relativistic approximation of a sum of ladder diagrams 

^In the approximation of annihilation at rest all the diagrams can be expressed as a linear combination of 
integrals with only two Feynman parameters, because in this case there are only two independent external 
momenta. 

® Except that we have to include the IR divergence of the W wave- function renormalization due to the 
photon exchange, which is cancelled by a real photon emission, see section 2.2.4. 

''The propagators and vertices of these diagrams give three powers of momentum in the numerator 
and eight powers in the denominator therefore the integration in four dimensions is convergent by power 
counting. In the case of diagram 5 the analytic integration on one parameter has been done using the 
PrincipalValue prescription, in order to discard the absorptive part, due to intermediate W^s being possibly 
on-shell. This part does not interfere with the tree diagram and thus would give a higher order contribution. 

The diagram containing the four vector boson vertex gives a vanishing contribution for the Wino anni- 
hilation at rest in a spin-singlet state. 

^Divided by the y/2, due to the difference in the normalization of the initial two-body state for neutralino 
and chargino pairs, see Eqns. (1.6) and (1.7). 



-8- 



having 1^-propagators as steps. Schematically: 

oo 

A = Y,^ox{W,tepr, (2.2) 

n=0 

where x means a non-relativistic convolution in momentum space and (VFstep)"" = ^step x 

• • • X Wstep- 

Since the radiative corrections of (VFstcp)" have been already treated, one has to com- 
pute the radiative corrections of Atree, that is replace 

Atree^ Atrcc + ^Di, (2.3) 
i 

where -D, are the various diagrams correcting ^treo but of course one must not include 
the contribution of non-relativistic approximation of the ladder diagrams that has been 
already taken into account. The contribution of the diagram 4 is then: 

Di = Atree X W^step + ^treeJ^C4. (2.4) 

2 

At a one- loop level Aree x Wstep = ^treef^;^, and therefore 

A r - n A 

(47rj^ 47r m\Y 

has to be considered as the genuine contribution to the radiative correction to ^trco- 

Diagram 5 is due to the exchange of Z or 7 between the final (see figure 2). This 
loop is IR divergent in the part in which there is the photon exchange. This is the only 
IR divergence in the radiative corrections of fxo; because the initial does not couple to 
the photon and there is no photon contribution to the wave-function renormalization, 
moreover in the other loops at least one of the x line is off-shell thus avoiding IR divergences. 

2.2.3 The total result for the radiative corrections due to the loops 

On figure 3 we present the results for the one-loop corrections coming from all the diagrams 
separately, in function of the DM mass being in range of 100 GeV - 3 TeV. One can see 
that the largest contributions come from diagrams 2,3 and 5, all of which contain photon 
exchange. Although the C5 is IR divergent, we get a finite result by giving (in all the 
numerical calculations) a small (with respect to the TeV scale) mass = 0.1 GeV to the 
photon. In reality it is of course massless, and indeed as we shall see the dependence on 
rriy will drop out in the final result. We will come back to this point in section 2.2.4, where 
we discuss the cancellation of the IR divergence by the inclusion of a real production. 

Actually, in order for the cancellation to be exact, we have also to take into account 
the IR divergent part of the virtual photon contribution to the W wave-function renor- 
malization, that is not included in the renormalization of g at the scale mw (see figure 
4). It gives a further contribution to C which is: s^41og(^). Including it in one- loop 
corrections gives finally 

Ci-ioop = 5^ + 44 log . (2.6) 
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Figure 3: The results for the one-loop correction to the amplitude of the x^X^ ~^ W^W' annihi- 
lation. The total correction is obtained by summing all those contributions and including the real 
production. The C5 contribution is made finite due to adding a small mass to the photon — 0.1 
GeV. In these result all the multiplicities of the diagrams were taken into account. 




Figure 4: The IR divergent diagram present in the gauge boson wave-function renormalization. 




From the computation we get that Ci^ioop < and also that it diverges for 0. 
The divergence is due to the graph 5 described above and to IR divergent part of the W 
wave-function renormalization. 

2.2.4 The radiative correction due to the real production and the cancellation 
of the IR divergences 

As we already mentioned, doing the computation of ao at the order 0{g^)^ we have also to 
add to (7o the production cross-section of W^W~Z and W^W~j^ which are of the same 
order (see figure 5). This has to be done on the level of the cross section. Hence, we start 
from a short review of the cross-section computation. If we call the amplitude A^, then 
the formula for the differential cross-section in our case reads: 

= (n wk^) E \M\WiP - E '=.). i^-n 

\ i / pol i 
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where Vr = 2v is the relative velocity, P = (2m, 0, 0, 0) and m is the mass of annihilating 
DM particle. The sum over polarizations gives: 

E'.< = -9,. + ^. (2.8) 

pol ^W,Z 



for massive gauge bosons and 



pol 



for the photon. 

In the annihilation into two particles with the same mass rrig the integration over the 
phase space gives: 



1 / TH 

' pol 

For the annihilation into three body final state, in the limit in which initial particles 
are in rest, the cross-section can be computed in a convenient parametrization with the 



use of Dalitz variables^, m?- = {ki + kj)'^: 



^^^3 = (n^ys^a(l H ~ E \M\^dmj^dml,. (2.11) 

pel 

The integration limits on these variables depend only on the masses and can be conveniently 
presented as 



4mi < < (2m - ma)^ {rr?2'^m.in < < {ml.^)max-, (2.12) 

with 

{ml^Uin = {E2 + E^f - (^^El-ml + , (2.13) 

{ml^)ma. = {E2 + E;f - I^^El -m\-^El- mij . (2.14) 

Here E2 = myijl and E-^ = {Am? — — mf2)/2mi2 are the energies of particles 2 and 3 
in the mi2 rest-frame. 

In order to add these contributions to the one-loop corrections we define the coefficients 



C^^ and C^^ as: 



®In actual numerical computations we follow a more direct approach by integrating over the final energies, 
which we check to be equivalent; this is numerically more convenient but the formulae are too long and we 
don't write them here. 
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t h 



b) 
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Figure 5: The diagrams for real production of tire gauge bosons and tlie production of the tb quark 
pair. 



Note the factor 2 in the definitions, which makes these coefficients to be the corrections to 
the amplitude coming from the real production. Using this we can write that the production 
cross-section provides a further correction: 

— loop+rp 

The full Ci 

-loop+rp coefficient should go to a finite constant for — >■ 0, which we find 
it is indeed the case. On the right plot of figure 6 we show the separate contributions from 
one-loop corrections and the real production to show that their sum is independent of m^. 

Three body production involving t quark. There is a further, though very small, 
contribution to the production cross-section at the order 0{g^): the processes involving 
the t quark in the final state 

x°x° ^ w-td, xV° ^ w-ts, ^ w-tb, 

^ W+id, ^ W+is, x°X° ^ w+ib 

These processes are due to the couplings — )• td, ts, tb and their conjugate. 

Notice that the other processes where in the final state there are a charged W and 
either a charged klj or a lighter charged qiQj pair must not be included, because they sum 
up to the total width of the charged W, and therefore are implicitly taken into account 
by unitarity when one takes the approximation of considering W as a stable particle. But 
the top is more massive than the W and therefore the charged qq pairs where one of the 
quarks is t are not included in in the total width of the W and have to be added to the 
correction. 

Since the square of the coupling — t- td is negligible with respect to the sum of the 
square of — )■ tb and — )■ ts (which all together add up to 5^), and the masses of b, s 
are negligible at our energy scale, by defining as before 

^W?^' ^ — — ' (^-^^^ 

we get final result for the total correction to the tree amplitude 

C^.ioop^rp+t = V + 4 f 4 log f — ^ + C'-A + clC'i + Ct . (2.18) 
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Figure 6: Left plot: the correction to the tree level x^x" ~^ W^W" amplitude coming from loop 
corrections (dashed red line), real production (chain blue) and the t quark (dotted violet). The full 
result is given by the solid black and sparse green lines (without and with the one-loop Sommcrfeld 
correction, respectively). Right plot: the dependence of the full result on the photon mass, for fixed 
m = 1 TeV. A complete cancellation of the IR divergent terms can be seen, and that the full result 
is independent of m^. 



In the numerical results, as we will see, the relative contribution of Cj is very small, 
due to the fact that it does not contain any large Log's, which are present in the case of 
the production of three gauge bosons. 

Note, that since the unphysical neutral Higgs is not coupled to W+ we can use the 
Feynman-'t Hooft gauge for the vector bosons forgetting the unphysical Higgs. As for the 
physical Higgs, its coupling to W+ is proportional to gmw and and therefore the (virtual 
or real) process involving it will be suppressed by a factor m^^ / rn? and we neglect them. 

2.2.5 The total result for the annihilation of x^X*^ 

We show the results for the full radiative corrections to the annihilation amplitude 

on the left panel of figure 6 (by now without the full Sommerfeld effect). One can see that, 
subtracting the one- loop Sommerfeld effect (that will be non-perturbatively treated), the 
total corrections (the solid black line) are significant, reaching over 15% for the m = 3 
TeV, but still in the perturbative regime. 

We also see that at a TeV scale the one-loop perturbative evaluation of the Sommerfeld 
effect is quite large and this is one of the reasons why the full non-perturbative treatment 
of this effect is needed. 

In writing a = 5 x ctq we are assuming that, while S is the non-perturbatively evaluated 
Sommerfeld effect, we can compute ctq by the standard Feynman diagrams of perturbation 
theory. However, for growing m/miy the contributions to (Jq of the diagrams 0{g^) becomes 
larger and larger, and the perturbative evaluation of ctq looses its meaning. Indeed, the 
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Figure 7: The diagrams for the one-loop corrections to x~^X ~^ W~^W annihilation. 

perturbative evaluation of the correction to do looks like to be border-line-reliable up to 
values of m of a few TeV. 

This fact is not surprising: when m and therefore the overall scale of the process gets 
large as compared to mw, the vector bosons resemble more and more to massless would-be 
gluons of an unbroken SU{2), like an SU{2) version of QCD. There occur large Log's of 
the ratio m/mw, and powers of them, which are not related to the UV divergences (and 
therefore cannot be included in a standard renormalization group treatment). Therefore, 
for higher values of m, one would need to borrow from QCD sophisticate techniques of 
re-summation of powers of large Log's or semi-empirical formulae. All that is beyond the 
scope of this work. 

2.3 The radiative corrections to x^X annihilation 

Due to the Sommerfeld effect the x^X~ annihilation gives a non-negligible contribution to 
the x^X^ annihilation process, which in fact can be of the same order as the direct process. 
Therefore, it is also important to compute the radiative correction to annihilation with 
X~^X~ ill the initial state. Because the computations are very similar to the x^X^ case, we 
don't discuss all the computations in detail, but rather stress the differences and present 
the final results. 

In this case, in the Feynman-'t Hooft gauge it occurs also the vertex of the charged 
unphysical Higgs with the vector bosons. However, in the same way as for the physical 
Higgs, its coupling is proportional to gmw and therefore the process involving it will be 
suppressed by a factor rn^f/m? and we neglect it. 

2.3.1 One-loop corrections to x^X ~^ W^W~ 

In the case of the annihilation of x'''X~ ) since they are charged, there are more diagrams 
to be computed, see figure 7, the technique is however exactly the same. 
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Figure 8: The diagrams the correction to the process x^X 
production of Z, 7. 
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Figure 9: The correction to the x~^X ~^ W~^W ampUtude (left plot) and dependence on the 
photon mass (right plot). The notation is the same as in figure 6. 



Note also that due to the difference in the normalizations of the initial states (Eqns. 
(1.6) and (1.7)), at the tree level A^^.^^^^^ = j.A'^S'^,^^+w- 

2.3.2 The radiative correction due to the real production 

Also in this case the computations goes in the same way, except that now the initial state 
particles are coupled to Z and 7, which gives the initial state Bremsstrahlung process 
(instead of internal one as in the x^X^ case). The diagrams to be computed are those on 
figure 8. 

2.3.3 The total result for the annihilation of x^X~ W^W~ 

On figure 9 we show the full radiative correction to the amplitude of the process x^X~ ~^ 
W~^W~ . When compared to the case of nuetralino annihilations, one immediately sees 
that although results are qualitatively similar, quantitatively are considerably smaller. In 
fact, the full one- loop result without including the one-loop Sommerfeld effect is within 
-10% range even up to 3 TeV. 




Figure 10: The diagrams for the one-loop corrections to x^X annihilation to neutral gauge 
bosons. 

For completeness, also in this case we show the IR cancellation and that our results 
are independent of m^. 

2.4 The one-loop corrections to x^X~ ~^ ZZ, Z7, 77 

The diagrams to be computed are given on figure 10. In this case there is no wave-function 
renormalization of the final states, because they do not couple to the photon and thus do 
not exhibit IR divergences. Moreover, in this case there are no IR divergences in the total 
one-loop corrections, since the fermion wave- function renormalization cancels precisely the 
IR divergence coming from the correction to the initial states (the bottom left diagram). 

There is also no three body production, since the emission of three (a mixture of 
Z and 7) is forbidden by the CP conservation: the initial state being spin singlet has an 
even CP, while both Z and 7 are CP-odd. ^'^ 

The results for the radiative correction to these processes are presented on figure 11. 
The corrections are very similar to each other, as could be expected from the fact that 
since m is much larger than m^, the differences in masses of the final states are not very 
important. On the other hand, the differences in couplings are taken into account in the 
tree level amplitudes for these processes (i.e. every of these three corrections is normalized 
to its own tree level amplitude). 

One can also see that the absolute value of these corrections is quite large, in fact 
considerably larger than for the annihilation into charged final states. This might look 
surprising, since there are less diagrams and none is IR divergent, but actually it can be 
easily understood by the fact that in this case there is no compensating effect of the real 
production. 

^°The processes involving two neutral gauge bosons one of them subsequently decaying into quark or 
lepton pairs are allowed, but similarly to what was said for the t quark production they are very suppressed 
and therefore we neglect them. 
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Figure 11: The full one-loop corrections to the annihilation of x^X~ i^ito ZZ,Z^,^^. These 
corrections come only from the loop diagrams, because there is no real production in this case (see 
the text). 



3 Final result for the Sommerfeld enhanced av 

Finally, we are ready to discuss the total annihilation cross sections of the x^x^ pair 
including both the radiative corrections discussed in the previous sections, as well as the 
full non-perturabative Sommerfeld enhancement. 

The results are presented on figure 12. Let us first concentrate on the left plot showing 
the x^x^ — )• W~^W~ (including the three body production) annihilation process. Our full 
results including the one-loop corrections and the Sommerfeld effect are given by the solid 
black line. There is a clear resonance visible, which is due to the creation of a loosely 
bound state of the incoming neutralino pair. The resonance occurs approximately when 
the Bohr radius of the X^X^ P^^ir matches the interaction range, i.e.: 

« . 3.1 

ma mw 

For this reason, the position of the resonance depends strongly of the value of the coupling 
and this is why it is so important to use it at an appropriate scale. In section 2.1 we have 
discussed that the proper value is the one at the electroweak scale. If one uses instead its 
runned value at the scale m (and do not include the radiative corrections) then one get the 
result plotted by the brown dotted line. That is, if one used the running of the couplings 
instead of doing full one-loop computation, one would get the resonance peak displaced in 
m from about 2.35 to about 2.5 TeV. The shift of this peak influences also the mass of the 
Wino for which the correct thermal relic density is obtained (see ref. [22]). 

On the same plot we show that including only the one-loop approximation of the 
Sommerfeld effect is a good approximation only up to about 200 GeV, beyond which it 
breaks down, mainly due to the presence of the resonance. However, it is of course still 
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Figure 12: Left plot: the total cross-section for tire anniiiilation of x°X° to W~^W~ (including 
the three body production). Our final results including both the one- loop corrections and the 
Sommcrfcld effect (SE) arc given by the solid black line. For comparison we plot the tree result 
(solid blue), tree level with the full SE (chain red), full one- loop level results but without non- 
perturbative SE (twin green) and the tree level with (dotted brown) and without SE (sparse blue) 
but with runned couplings at the scale m. Right plot: the cross section for the annihilation to 
ZZ, Z"f,^^. The full one-loop results with the SE included are given and for comparison the 
leading order ones. 



more accurate than just using the tree level value. What might seem surprising is that 
using the tree level formula but with a running coupling constant at a scale m (blue sparse 
line) is even a worse approximation than simply taking the standard tree level cross-section. 
This comes from the fact that running of the couplings captures only the UV effects of 
re-summation of large Log's (which gives a negative contribution), while in our setup the 
dominant correction to the annihilation amplitude is the Sommerfeld one, which is positive 
and (at a one-loop approximation) proportional to m/mw- 

We also do a comparison with the results using the full Sommerfeld corrections, but 
this time applied only to the tree level annihilation amplitudes (red chained line). The 
largest difference occurs before the resonance, for the masses of about 1 TeV, where the 
inclusion of one-loop contributions makes the full cross-section smaller by as much as about 
30%. In the resonance, although due to the logarithmic scale from the plot it seems that 
including the loop corrections do not change the result significantly, actually the radiative 
corrections make the value of av in the peak is smaller by about 22%. 

On the right plot of figure 12 we present the results of the annihilation cross section 
to ZZ, Z7 and 77. Again a clear resonance is visible, for the same reason as before (in 
fact, by construction Sommerfeld effect is independent on the final states). However, the 
absolute value of the cross-section into 77 is about two orders of magnitude smaller than 
into W~^W~ . Since the annihilation of x^X^ ™to neutral gauge bosons cannot occur at the 
tree level, the Leading Order (LO) for those processes is computed by including only the 
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one-loop contributions (from which the one-loop Sommerfeld is dominant for large masses, 
which makes the LO proportional to rather than m~^, what explains the shape of 

the LO results). 

Looking at the low mass region, one can see that the full cross-sections are smaller 
than the LO ones. There are two reasons for that, both relaying on the fact, that, as 
we discussed, the one-loop corrections are negative. This means that both annihilation 
amplitudes, to be combined with the Sommerfeld enhancement factors, are decreased with 
respect to their tree level values and, since the annihilation x^X^ to neutral gauge bosons 
do not occur at tree level, the two terms from which the Sommerfeld enhanced amplitude is 
constructed (see eq. (1.14)), have opposite signs. Thus, not only the one-loop contribution 
decrease the x'^X~ annihilation, but also the Sommerfeld effect is suppressing the cross- 
section when So ^ s±, i.e. precisely in the low m region. 

In the Appendix A we provide a fit to our results, listing in a Table 1 the numerical 
values for the parameters. This fit could be of use for possibly incorporating our results in 
further studies of the dark matter signals. 

Finally, let us make explicit which perturbation theory orders are included in our 
computation. The resulting cross-section a2 can be written as (see eq. (2.10)): 



a2V = 




where A^^ and are the amplitudes of annihilation into two gauge bosons for initial 

X°x*^ and x^X ) respectively. The subscript refers to the order of perturbation theory, i.e. 
A2 is the 0{g^) part and ^4 the 0{g^) one. The Sommerfeld factors sq and s± (see eq. 
(1.13)) are evaluated non-perturbatively. 

One can see that we do not include any terms of the order 0{g^), say |^4p, or higher. 
They are subdominant at the TeV scale, but may play a role in the low mass region 
(especially for the annihilation to ZZ, 77 and Z-y, where the tree level is not present). 
In order to incorporate them consistently with the Sommerfeld effect one would need to 
include also additional corrections, e.g. two-loop contribution to the x^X~ annihilation, 
which is beyond the scope of our paper. 

4 Conclusions 

In this paper we have computed the full annihilation cross-section including one-loop and 
the Sommerfeld electroweak corrections for the annihilation of the non-relativistic Majorana 
fermion being in the adjoint representation of the weak SU (2). This scenario can be realized 
in number of physical situations, from which the most important is the neutralino in the 
MSSM in the limit in which it is purely Wino. The most interesting range of masses in this 
case is about a few TeV, because then the Wino can account for the thermal dark matter 
with the relic density in agreement with the observations. However, for such large values of 
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the mass the perturbative computations are not sufficient. To give a correct annihilation 
cross-section in this case one needs to sum over all orders the ladder diagrams which give 
the Sommerfeld enhancement. 

In our work we included both these effects and discussed how to incorporate them 
simultaneously. In order to do so, one needs to compute both the x^x" ^-iid X^X' annihi- 
lation amplitudes at the same order, multiply them by corresponding Sommerfeld factors 
to get the corrected amplitude and finally take its modulus squared and integrate over the 
phase space. When computing the one-loop corrections one also has to be careful and do 
not include the one-loop version of the Sommerfeld correction, since it is included in the 
full non-perturbative treatment. 

The results we obtained introduce a relatively large one-loop corrections both to the 
X^lxP — )■ W^W~ and x^X~ ~^ W~^W~, ZZ, Z7, 77 processes. For a neutralino mass being 
m = 3 TeV the corrections to av reach more than -30% (excluding the contributions of 
the Sommerfeld effect) for the neutralino annihilation, -20% for the x~^X~ ~^ W~^W~ and 
more than -70% for charginos annihilating into neutral gauge bosons. As we discussed 
in section 2.2.5, this is due to the occurrence of powers of \ogm/m\Y , not related to UV 
divergences, that make the gauge theory to resemble a confining unbroken SU{2). 

On the other hand the Sommerfeld effect introduces a positive contribution, which 
already for masses of about 200 GeV overwhelms the other corrections and become far 
dominant when we enlarge the mass to the TeV scale. On top of that, the Sommerfeld 
factors exhibit a resonance due to forming a loosely bound state of the incoming particles. 
This has been already discussed in detail in literature (see e.g. [18-20, 31]) and we obtain 
qualitatively the same results. However, since in this work we included rather large one- 
loop contributions, the precise values of the cross-section are slightly smaller, also in the 
resonance. Our final result for x^X^ ~^ W~^W~ cross-section before the resonance is up to 
about 30% smaller than the one computed from the tree level cross-section multiplied by 
the Sommerfeld factor, whereas in the resonance it is smaller by about 22%. 

We also discuss the issue of the values of the couplings which we have to use, and 
argue that those should be the ones at the electroweak scale. This introduces a very mild 
modification to the one-loop corrections, but leads to a visible shift in the mass of the 
neutralino for which the resonance in the annihilation cross-section occurs (as well as in 
the relic density computations [21-23, 31]). Our results show that the resonance takes 
place for ni ~ 2.35 TeV whereas if we used the couplings runned up to the scale m then 
this would shift to about 2.5 TeV. 

The possible phenomenological applications of these computations are mainly concen- 
trated on the modification of the cosmic ray spectra and therefore the predictions for the 
indirect detection. There inclusion of one-loop contributions introduce two effects. First, 
the overall cross-section being lowered by the one-loop corrections, both in the annihila- 
tions into W~^W~ and ZZ, as well as Zj and 77. The two former processes contribute to 
the diffuse gamma ray spectra, while the two latter give a gamma line. Since our results 
show that the one-loop corrections are larger for the neutral final states than the charged 
ones, it suggests that the gamma lines are even fainter with comparison to the diffuse 
background. On the other hand, the second consequence of the full one-loop computation 



-20- 



is that now the total cross-section incorporates also the real production, which may alter 
the final spectra significantly (see e.g. [7, 8, 10, 11]). We leave the full discussion of this 
point to a future work [32]. 
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A Fit to the full Sommerfeld enhanced av 

In order to make our results easy to use in further computations, we also provide a fit 
of the full cross-sections (comprising the three body production) including the one-loop 
corrections and the full Sommerfeld effect. All of them are fitted with the same functional 
form: 



The values for the coefficients providing the best fit are given in Table 1, for the various 
annihilation channels. Both /<.>(?ti), with m and expressed in GeV, give the result 
in units of ctd? /s. The difference between these fits and the numerical results are within 
1% range (see figure 13). 




(A.l) 



where 




(A.2) 
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-6.24676 X 10-20 


-8.02201 X 10-21 


9.73209 X 10-21 


bo 


bi 


b2 


bs 




2.67711 X 10-^2 


1.35313 X 10-1^ 


-7.00791 X 10-19 


2.49388 X 10-22 




Co 


Cl 


C2 


C3 




-5.26864 X 10^3 


1.44622 X 10" 


-9.46014 X lO'^ 


1.81301 X 10"^ 




x°x° ^ zz 


ao 


ai 


a2 


as 


a4 


-1.15577 X 10^27 


5.20277 X 10-27 


1.86739 X 10-27 


-2.80783 X 10-27 


7.87796 X 10-28 


bo 


bi 


b2 


bs 




3.3929 X 10"^^ 


-7.98612 X 10-19 


-1.2351 X 10-22 


4.67818 X 10-26 




Co 


Cl 


C2 


cs 




8.52321 X 10" 


-8.39761 X 10*^ 


2.51261 X 10^ 


-2.05257 X lOi 






ao 


ai 


a2 


as 


a4 


8.91573 X 10^28 


1.11734 X 10-27 


2.27249 X 10-27 


-1.99817 X 10-27 


5.12998 X 10-28 


bo 


bi 


b2 


bs 




-1.80836 X 10^1^ 


3.37165 X 10-1^ 


-7.40609 X 10-22 


7.48779 X 10-26 




Co 


Cl 


C2 


Cs 




-7.32662 X 10" 


1.74234 X 10^ 


-1.0862 X 10^ 


2.03233 X 102 




X°X° 77 


ao 


ai 


a2 


as 


a4 


3.64069 X 10^2^ 


-1.73682 X 10-28 


5.95796 X 10-28 


-3.91584 X 10-28 


9.08957 X 10-29 


bo 


bi 


b2 


bs 




4.00871 X 10-12 


-2.48789 X 10-15 


5.60019 X 10-19 


-5.44691 X 10-23 




Co 


Cl 


C2 


Cs 




8.04594 X 10^^ 


-9.37744 X 10i2 


3.60859 X 109 


-4.57474 X 10^ 





Table 1: The obtained values of the parameters giving the best fit for the cross-sections crvf^t for 
all y^x^ annihilation channels. 
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Figure 13: The comparison between tlie numerical results (points) and the fit (solid lines) for 
the annihilation cross-section of x^x" to W^W~ (top) and neutral gauge bosons (bottom). The 
accuracy of the fits is within 1%. 
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